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We study the magnetic susceptibility of ID quantum XY model, and show that when the temperature ap- 
proaches zero, the magnetic susceptibility exhibits the finite-temperature scaling behavior. This scaling behavior 
of the magnetic susceptibility in ID quantum XY model, due to the quantum-classical mapping, can be easily 
experimentally tested. Furthermore, the universality in the critical properties of the magnetic susceptibility in 
quantum XY model is verified. Our study also reveals the close relation between the magnetic susceptibility 
and the geometric phase in some spin systems, where the quantum phase transitions are driven by an external 
magnetic field. 
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Introduction: Quantum Phase Transitions (QPTs), which occur at absolute zero and are driven by zero-point quantum fluc- 
tuations, are one of the most fascinating aspects of many-body systems. QPTs and related quantum critical phenomena have 
been a topic of tremendous interest in condensed matter physics and have been extensively studied in the past decade |[T). In 
recent work, quantum criticality has been characterized by using the methods and notions borrowed from quantum information 
science, such as the concurrence |2|, the entanglement entropy [3 1, geometric phase (GP) [4], Loschmidt Echo Q, and quantum 
fidelity [ 6 1 in the place of traditional criteria, such as specific heat or magnetic susceptibility (MS). Most of these studies focus 
on the zero-temperature properties of the critical systems. In recent years, the finite-temperature properties of QPTs Q0> such 
as, thermal entanglement [9 | have begun to attract more attention. This is because, firstly, all experiments are confined to finite 
temperature. Thus, to experimentally verify the theoretical results, knowing only the zero-temperature properties of the quantum 
system is not sufficient. Secondly, though genuine QPTs occur only at absolute zero, quantum criticality has profound influence 
on system properties up to a surprisingly high temperature [7 1. Interesting phenomena at finite temperature related to QPTs have 
been experimentally observed in various systems, such as the heavy fermion system and the BEC ifTOll . 

On the other hand, it has been shown that a QPT in d space dimensions is related to a classical transition in d + z space 
dimensions [T] QT), where z is the dynamical critical exponent. Under this quantum-classical mapping, the temperature T 
of quantum system maps onto an imaginary time direction: t = —it/h 6 (0, l/fc^J 1 ), where r and t are imaginary and 
real time ATI . Accordingly accessing the QPT by reducing the temperature amounts to increasing the size of imaginary time 
dimension towards infinity, and leads to a divergence of the spatial correlation length £. This one-to-one mapping motivates 
us to study the finite-temperature properties of QPTs through its higher dimensional classical counterpart. Studies of these 
QPTs and the quantum-classical mapping rely heavily on the exactly solvable models. One of the most common examples is 
the one -dimensional quantum Transverse Ising Model (ID TIM) lfj"2l . which exhibits a second-order QPT at the critical point 
A c = 1, and its classical counterpart - the two-dimensional classical Ising model [13], which exhibits a second-order thermal 
phase transitions at the Curie point. 

Historically, scaling has played a central role in the study of classical criticality. It is well known that the 2D classical Ising 
model obeys finite-size scaling behavior fl4l . A straightforward idea is to study T ^ scaling laws of ID TIM. In Refs. |9| 
and 1151 . the authors use Gmneisen Parameter and concurrence to characterize finite-temperature properties induced by QPT at 
zero temperature. In this paper, instead we will use a classical macroscopic thermodynamic obsevable - the MS - to study the 
finite-temperature properties of the generalized ID TIM - the quantum XY chain. The MS has the advantage of being easily 
experimentally accessible and has been used as a witness of macroscopic quantum entanglement ifTrjlfTTll . We will show how 
the finite-temperature scaling is manifested when the temperature approaches zero, in analogy with finite-size scaling in the 
imaginary time direction of the 2D classical Ising model. We will also verify the universality in the properties of the MS in 
quantum XY chain. Finally we will elucidate the close relation between the MS and another well studied observable - the GP 



Magnetic susceptibility of quantum XY chain at finite temperature: The Hamiltonian of quantum XY chain can be written as 



where N is the number of spins in the chain; J is the coupling strength (for simplicity we choose J = 1 hereafter); A is 
external magnetic field, and 7 describes the anisotropy of the system; erf, a = x,y,z are the Pauli matrix on the ith site of 
the chain. After a standard procedure fl"2l . this Hamiltonian can be diagonalized as H(j, A) = J^k ^^kiwlvk ~ 1/2), where 
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rjk is the Fermionic annihilation operator of the fc-th mode quasi particle; = \J{X — cos k) 2 + j 2 sin 2 k are one half of the 
excitation energy for modes k = 2n(i — 0.5) /JV, i = 1, 2, • • ■ , JV/2. The partition function of the system can be obtained as 
Z = Yik (e _/3Afc + e^ Ak ) = Ylk 2 cosh (/3Afe), where (3=1/ fesT is the inverse temperature and k B is the Boltzmann constant. 
Accordingly, the free energy per spin of the system can be calculated as J 7 = —ksTln Z/N = — k B T m P cosh (/3A&)] /JV. 
In the thermodynamic limit, JV — > oo, we use an integral to replace the sum and obtain the exact expression of the free energy 
per spin at temperature T 1 12 1 

i r 

F = ~k B Thi2- k B T x - dfcln[cosh(/3A fc )] . (2) 

f Jo 

The magnetization per spin along the direction of the external magnetic field A at temperature T can be obtained 

M Z {T) =~ = - f tanh(/3A fe ) X ~™ sk dk, (3) 

OA 7T J A k 

and then the MS along z direction Xz = —d 2 F/dX 2 as a function of the temperature T and the magnetic field A of the system 
can also be obtained 
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We plot the MS Xz of ID TIM (7 = 1) as a function of external magnetic field A and the temperature T in Fig. 1. Clearly it 
can be seen that the logarithmic divergence of the MS at zero temperature indicates the second-order QPT at the QCP A c = 1. 
We would like to point it out that at zero temperature, the magnetization is reduced to M Z (T = 0) = J Q (A — cos k) /(%Ak)dk. 

For the convenience of later study, we introduce another observable - the GP, which is a fundamental concept in quantum 
mechanics [21]. To obtain a geometric phase, we rotate the Hamiltonian ([I} around the z axis for an angle </>. The effective 
Hamiltonian after the rotation is 

iV 

= U^HUl U+ = J] e^/ 2 . (5) 
i=i 

The periodicity of the Hamiltonian in <j> is tt. After we rotate the Hamiltonian back to its initial form (</> = 71"), the GP of the 
ground state accumulated by varying the angle <f> from to it is given by 

P g = -i^£ ((GS\UI) ^{U^\GS))d^ (6) 

which is an extra phase in addition to the usual dynamic phase. From Refs (4J [T8] [19] [20] we know that the ground-state GP 
studied there can be expressed as 

Pg = 7T + r A ~° OSfc rffc = 7T + nM z (T = 0). (7) 

Jo A k 

Hence, the derivative d(3 g /d\ of the ground-state GP over the external field is ir times of the zero-temperature MS Xz — 
dM z (T = 0)/d\. We can understand this relation in the following way: the ground-state GP studied in Refs. ]4l [181 [T9"l 
l20l is a function of the derivative of the ground state energy with respect to the external magnetic field |[T8l [T9l l20l . and at 
zero temperature, the free energy is equal to the ground state energy. Thus, at zero temperature, the GP is a function of the 
magnetization. 

As is well known, at zero temperature, the MS of ID TIM shows logarithmic singularity at the QCP and exhibits finite-size 
scaling behavior in the proximity of the QPT point A c = 1. Thus, it is not surprising that the GP exhibits singularity and 
finite-size scaling behavior near the QCP 1 18 1. Instead of studying the finite-size scaling of the GP (MS at zero temperature), in 
this letter, we will study finite-temperature scaling of the quantum XY chain. We will see that when the temperature approaches 
zero, in analogy with the imaginary time direction approaching the infinity in the finite-size scaling, the MS obeys T / scaling 
behavior in the proximity of the QPT. 

Scaling of the magnetic susceptibility of the quantum XY chain: In order to further understand the relation between the ID 
TIM and 2D classical Ising model, we investigate the finite-temperature scaling behavior of the MS by the finite-size scaling 
ansatz [22 1 . For simplicity, we first look at ID TIM (7 = 1), and we will discuss the properties of the family of 7 7^ 1 later. 
The MSs as a function of the external magnetic field A at different temperatures T (including zero temperature) are presented 
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FIG. 1: MS Xz °f ID TIM model as a function of external magnetic field A and temperature T. It can be seen that the MS at zero temperature 
show logarithmic divergence at the QCP A c = 1. At nonzero temperature the MS is analytical. This agrees with the known result that ID TIM 
model does not exhibits thermal phase transition at nonzero temerature. 




FIG. 2: (color online), (above) The MS for the ID TIM (7 = 1) as a function of the controlling parameter A. The curve corresponds to different 
temperatures ksT — OJ, 0.02J, 0.06J, 0.21J, 0.5J, and 1.01J. With the decrease of the temperature, the maximum gets pronounced, and 
(below) the pseudopoint A m changes and tends as T 1706 towards the QCP A c = 1. 



in Fig. 2. At zero temperature the MS shows a singularity at A c = 1, but at nonzero temperature, there are no real divergence 
of Xz- Nevertheless, there are clear anomalies at low temperature, and the height of which increases with the decrease of the 
temperature. This can be regarded as the precursors of the QPT. What is more, the position A m of the maximum susceptibility 
(pseudocritical point) [22| changes and tends as T 1 - 704 towards the QCP and clearly approaches A c when T — > (see Fig. 2b). 
Meanwhile, the maximum value Xz\x m of the MS diverges logarithmically with the decrease of the temeprature 

Xz\\ m ~ Ki InT + const. (8) 

Our numerical results (see Fig. 3a) give k\ = 0.320. On the other hand, when T = 0, from Ref. 11231 we know that the MS in 
the proximity of the QCP exhibits logarithmic singularity 

Xz ~ K 2 In |A — A c | + const. (9) 

Our numericals in Fig. 3b give the result K2 ~ 0.317, while the exact result [23| gives K2 — l/n ~ 0.3183. We would like to 
point it out that the coefficient K2 here is the same as that in Ref. fl8l . where the author gives K2 ~ 0.3123 and our numerical 
result is closer to the exact result k 2 = 1/tt. According to the scaling ansatz in the logarithmic singularities, the ratio |k2/«i| 
gives the critical exponent v that governs the divergence of the correlation length £ ~ |A — \ c \~ v . In our case, v 1.009 ~ 1 is 
obtained in the numerical calculation for the ID TIM, which agrees well with the known result about ID TIM |[T2l . Furthermore, 
by proper scaling and taking into account the distance of the maximum Xz from the QCP, it is possible to make all the data for 
the value of F = 1 — exp [xz(A) — X«U m ] as a function of (A — A m )/T for different temperatures T to collapse onto a single 
curve (see Fig. 4). This figure contains the data for temperatures ranging from kgT = e~ 3 J, e _4 J, e _5 J, e _5 5 J. These 
results demonstrate that the MS does obey the scaling behavior as the temperature decrease to zero, in analogy to the lattice size 
approaching the infinity in the finite-size scaling cases. 
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FIG. 3: (color online), (above) The maximum value of the MS at the pseudocritical point A m of the ID quantum XY chain as a function 
of temperature T. The slope of the line is 0.317 (0.394) for 7=1 (blue) (7 = 0.8 (red)), (below) The MS at zero temperature diverges 
logarithmically in the proximity of the QCP A c = 1. The slope of the line is 0.320 (0.401) for 7 = 1 (blue) (7 = 0.8 (red)). The ratio of the 
two slopes (below and above) for a fixed parameter 7 is equal to the critical exponent v. Here v ~ 1.009 {v m 1.017) for 7 = 1 (7 = 0.8) is 
obtained. The numerical results agree with the scaling ansatz and the universality of the XY model. 



In the following we will study the universality of the critical behavior of the MS. It is well known that the anisotropic XY 
chain (7 9 (0, 1]) belongs to the ID TIM universality, while isotropic XY chain (7 = 0) belongs to the XX universality. For the 
ID TIM universality, v = 1, while for the XX universality, v = 1/2. We will show that the finite-temperature scaling behavior 
of Xz also manifests the universality principle - the critical properties depends only on the dimensionality of the system and the 
broken symmetry in ordered phase. To verify the universality principle of the XY model, we consider the case for 7 7^ 1. The 
asymptotic behavior is also described by Eqs. ([8]) and d9j. From Fig. 3 we see that for 7 = 0.8 numerical simulation gives 
Ki Rj 0.394 and K2 ~ 0.401, while the exact result ll23l should be «2 = (771-) -1 w 0.398. As a result the critical exponent for 
7 = 0.8 is v = |k2/ki| ~ 1.017, very close to the exact value v = 1. Moreover we also verify that by proper scaling, all data 
for different temperatures T but a specific 7 will collapse onto the same curve. The data for 7 = 0.8 are shown in Fig. 4. What 
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FIG. 4: (color online). The value of F = 1 — exp [\z (A) — \z | A m ] as a function of (A — A m ) /T for different temperatures (dots of different 
colors) ksT = e _3 J, e _4 J, e -5 J, and e -5 ' 5 J. For fixed 7 (here we choose 7 = 1 and 7 = 0.8), all data collapse on a single curve, which 
agrees with the finite-size scaling behavior. The critical exponent v = 1 can be obtained from this figure. 

is more, through a similar analysis to that in Ref. [ 18 1, we can directly extract the finite-temperature scaling behavior of the XX 
(7 = 0) universality class. It can be found that, at zero temperature T = 0, for the XX universality, the magnetization can be 
written in the following compact form. 

1_ |arccosA,(0< A< 1) 

1,(A>1) ' {V> 

Accordingly, the critical exponent v = 1/2 and z — 2 can be extracted from the MS \ z = \/2{l — (A — ► 1~) fl8l . which 

is different from the TIM universality (v — 1 and z — 1). When we change the anisotropy 7 from 1 to 0, we find the range of the 
validity of the quantum scaling ansatz in A (Eq. (|9j) shrink gradually. The leading term of the MS crossover from -i- In (1 - A) 

to 1/2(1 — A)~2. Hence, when < 7 < 1, the scaling belongs to XY universality, while when 7 = 0, the scaling belongs to 
XX universality. Finally, we also would like to point it out that our numerical result shows that the scaling behavior of Eq. ([8]) 
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can persist up to a temperature ksT w 0.5J. This result agrees well with that of Ref. Q. In addition, the crossover line in the 
region < A < 1 given by the MS \z is roughly T c ~ |A — X c \ uz , which agrees well with the result obtained in the analysis 
elsewhere 0. Hence the boundary of quantum critical scaling region can be confirmed by the behavior of the MS \z- 

Magnetic susceptibility and geometric phase: As we have mentioned before, the derivative of the ground-state GP discussed 
in Refs. Bl [181 [191 1201 is equal to tt times of the MS, and the finite-size scaling of the GP [18] actually represents the finite-size 
scaling of the MS. Based on these studies, we would like to further study the relation between the thermal-state GP and the 
MS at a finite temperature. Similar to the definition of the ground-state GP in Refs. (4] [T8] [T9j |20), we define the thermal- 
state GP in the following way: four eigenstates of the modes (k, —k) of (see Refs. JH [T8l QjO |20)) can be expressed 
as 100)* = co S (9 k /2)\0) k \0)_ k +i^sm{9 k /2)\l} k \l}_ k , \ll) k = ie^sm (9 k /2) |Q) fc |0)_ fc + cos (6 k /2) \l) h \l)_ k , 
\0l) k = |0> fe |l)_ fc , and |10) fc = |l) fc |0)_ fe with the angle 6 k defined by 9 k = arctan[— sin fc/(cos k - A)]. The GP of the 
thermal state at temperature T accumulated by varying the angle </> from to tt is described by 

_o • N/2 r Ft 

^-tEE/^h4^* (id 

fe=l n J V 

where \n) k = \00) k , |01) fe , \W) k , and |H) fc - After a straightforward calculation, we obtain the same relation between the 
magnetization and the GP as that of zero temperature /3t = 7r + Jq(^ — cos k)/A k tanh (0A k )dk = 7r[l + M Z (T)] (|3j. Thus 
we prove that at both zero temperature and nonzero temperature, the GP of the quantum XY chain is a linear function of the 
magnetization and the derivative of the GP is proportional to the MS. The discussions of the finite-temperature scaling of the 
MS in this letter can be alternatively regarded as the finite-temperature scaling of the GP in the proximity of the QPT point. 
Finally the close relation between the GP and the MS does not confined to ID quantum XY chain. In Ref. Il24l the ground- 
state GP of the Dicke model and its relation to quantum criticality are studied. We would like to point it out that, similar to the 
discussions about the ID XY chain, the ground-state GP of the Dicke model is a linear function of the ground state magnetization 
j3 g = 7r(l + (S x ) /N), where (S x ) /N is the magnetization M x per spin in Dicke model. Hence the derivative of ground-state GP 
of the Dicke model is also equal to tt times of the MS. Besides the above two examples, it can be proved that for any QPTs driven 
by an external magnetic field, such as the Lipkin-Meshkov-Glick model 1 19 1 and ID XXZ model |25 1, the relation between the 
GP and the magnetization still holds true. The proof is give as follows. For those QPTs driven by an external magnetic field, 
we apply a 7r-rotation along the z axis for every spin = Ylj=i e l ^" Tj ^ 2 to obtain the GP. The GP of the ground state can be 
expressed as (|6|l 

g = -ijj J* ({GS\UI) (i^^(U^\GS))d^=^^2{GS\ a l\GS)=nM z (T = 0), (12) 

where \GS) is the ground state of the Hamiltonian (JT|) before the rotation, and \GS) is the instantaneous ground state after 

the rotation for an angle <f>. We know that J2j (GS\ ^f- \GS) is the definition of the GP of the ground state. Thus we prove that 
the GP obtained by applying a rotation around the z axis to each spin, is proportional to the magnetization along z axis (similarly 
if we rotate along x axis, the GP will beproportional to the magnetization along x axis). We also would like to point out that the 
GP in Eq. ( 12 1 differs from that in Eq. ill by a constant tt. This is because the ground state of has an uncertainty of the 



global phase. When we choose a proper global phase, we can eliminate the difference between Eq. |7]i and Eq. ( 12 1. When we 



study the scaling of ^§S-, the difference does not affect. In addition, we can generalize the above discussions to eigenstates other 
than the ground state. We find the same proportional factor between the GP and the magnetization for all eigenstates. Thus the 
relation between the GP and the magnetization can be straightforwardly generalized to a thermal state at finite temperature. The 
above XY model is a good example. 

In summary, we study the finite-temperature scaling of the MS of the quantum XY chain. All key features of the quantum 
criticality, such as scaling, critical exponent, the universality, etc. are presented in the MS of the XY spin chain. Though the 
nature of the QPT and the T ^ scaling is purely quantum mechanical, the classical macroscopic thermodynamic observable 
MS, which can be easily accessed experimentally, can be used to witness and characterize the quantum features of the system 
|[T6l[T7l . Our studies shed light on the mechanism of bring quantum criticality up to a finite temperature, and opens the possibility 
of observing the footprint of quantum criticality experimentally. We also would like to point it out that the results obtained in 
this paper does not depend on the model and thermodynamic observable used here and can be generalized to other QPT models 
with the only change of MS to a "controlling parameter-dependent susceptibility". For example, in a QPT driven by the pressure 
instead of the external magnetic field, the observable \v = —d 2 F/dp 2 is expected to exhibit the finite-temperature scaling 
behavior, and the critical exponent can be extracted through a similar analysis. Finally, our study establishes the connection 
between the MS and the GP at both zero temperature and nonzero temperature in a family of spin systems, where the QPTs are 
driven by an external magnetic field. 
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